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We show that in a sufficiently strong magnetic field the QCD vacuum may undergo a transition 
to a new phase where charged p mesons are condensed. In this phase the vacuum behaves as an 
anisotropic inhomogeneous superconductor which supports superconductivity along the axis of the 
magnetic field. In the directions transverse to the magnetic field the superconductivity is absent. 
The magnetic-field-induced anisotropic superconductivity - which is realized in the cold vacuum, 
i.e. at zero temperature and density - is a consequence of a nonminimal coupling of the p mesons 
to the electromagnetic field. The onset of the superconductivity of the charged p mesons should 
also induce an inhomogeneous superfluidity of the neutral p° mesons. We also argue that due to 
simple kinematical reasons a strong enough magnetic field makes the lifetime of the p mesons longer 
by closing the main channels of the strong decays of the p mesons into charged pions. 

PACS numbers: 12.38.-t, 13.40.-f, 25.75.-q 



I. INTRODUCTION 

Properties of QCD matter subjected to very strong 
magnetic fields have recently attracted increasing interest 
of the community. The interest is motivated by the pos- 
sibility to create strong magnetic fields in the heavy-ion 
collisions at RHIC and LHC. The strength of the mag- 
netic field is estimated to be of the hadronic scale [JJ, 
eB ~ (1 . . . 15) m 2 , or even higher (here « 140 MeV 
is the pion mass). The duration of the magnetic field 
"flashes" is expected to be rather short (a few fm/c). 

Both analytical calculations [2H1] and lattice simula- 
tions [5 indicate that the QCD phase diagram is af- 
fected by the strong magnetic field. In particular, the 
external magnetic field splits the chiral and deconfine- 
ment transitions [31 3] . In a constant magnetic field of 
the typical LHC magnitude, eB ~ 15m 2 . [JJ, the split- 
ting between the critical temperatures of these transi- 
tions reaches 10 MeV |3j . 

In the quark-gluon plasma the strong magnetic field 
may also lead to the chiral magnetic effect 6] . This effect 
generates an electric current of quarks along the magnetic 
field axis provided the densities of left- and right-handed 
quarks are not equal. 

In the cold matter the external magnetic field may cre- 
ate spatially inhomogeneous structures which are made 
of quark condensates [7j. 

A recent lattice simulation has revealed that in the 
cold confinement phase the external magnetic field in- 
duces nonzero electric conductivity along the direction 
of the field, thus transforming the QCD vacuum from an 
insulator into an anisotropic conductor [8]. In our paper 
we argue that there is a chance that a stronger magnetic 
field may be able to make the QCD vacuum unstable 
towards creation of a superconducting state. We would 
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like to stress that we discuss here the electromagnetic su- 
perconductivity which should be distinguished from the 
color superconductivity in the dense matter [SJ. We dis- 
cuss a superconducting state which may presumably be 
formed in the cold vacuum, i.e. at zero temperature and 
density. 

Basically, we follow the works of Ambj0rn, Nielsen and 
Olesen on two subjects: (i) on the condensate of color 
magnetic flux tubes ("spaghetti states") [10] created by 
an unstable gluonic mode in the QCD vacuum [11]; and 
(ii) on the condensation of the VF-bosons in the stan- 
dard electroweak model due to sufficiently strong exter- 
nal magnetic field [H [13]. The key idea Refs. [TMS] 
is that the vacuum of charged vector particles is unsta- 
ble in the background of a sufficiently strong magnetic 
field provided these particles have anomalously large gy- 
romagnetic ratio g = 2. The large value of g guarantees 
that the magnetic moment of such particles is too large 
to withstand a spontaneous condensation at sufficiently 
strong external magnetic fields. 

As we have mentioned, there are at least two examples 
of such instabilities. A strong enough chromomagnetic 
field leads to the instability of the gluonic QCD vacuum 
since the gluon is the vector particle with the (color) gy- 
romagnetic ratio g = 2 [11] . As a result of the instability, 
a spaghetti of the chromomagnetic flux tubes is formed. 
These flux tubes tend to arrange themselves into a lattice 
structure similar to the Abrikosov lattice which is realized 
in a mixed state of a type-II superconductor subjected to 
a near-critical external magnetic field [TO] . 

The second example is suggested to be realized in the 
standard electroweak model. The gyromagnetic ratio of 
the W boson is also large, g = 2, so that in the strong 
magnetic field the vacuum of the electroweak theory is 
unstable towards formation of the condensate of the W 
bosons. The W condensate is accompanied by a similar 
lattice vortex state [12j [13] . Note that in the second 
example the external field is the electromagnetic field and 
not the color (gluon) one. 
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Our work is based on the fact that the p meson is the 
charged vector particle with the gyromagnetic ratio g = 2 
so that this particle may condense in a background of 
strong enough magnetic field. It important to stress that 
in all discussed cases of the spontaneous condensation 
- we mentioned the gluons in QCD [IT], the W bosons 
in the electroweak theory [T2] [T3] and the p mesons in 
QCD (this article) - the condensation takes place in the 
vacuum at zero temperature (as opposed to dense and/or 
hot environment). 

The structure of the paper is as follows. In Section |TT| 
we outline the basic idea of the p-meson condensation. 
In the same section we argue that the p mesons are (at 
least, partially) stabilized by the strong magnetic field 
background. This is an important property which should 
make the p condensate "intrinsically" stable against de- 
cays of the p mesons (the p mesons have a very short 
lifetime in the absence of the external fields). In Sec- 
tion III we describe the quantum electrodynamics of the 
p mesons. Section |IV| is devoted to a short overview 
of basic features of the Ginzburg-Landau model of the 
superconductivity (homogeneity, isotropy, effects of the 
magnetic field, the Abrikosov vortices, the Meissner ef- 
fect, the London equations). In Section [V] we discuss the 
same features in the superconducting state of condensed 
p mesons in QCD and find a few similarities and many 
surprising dissimilarities with the ordinary superconduc- 
tivity. The last Section is devoted to our conclusions. 



II. p MESONS IN STRONG MAGNETIC FIELD: 
CONDENSATION AND LONGER LIFE 



A. Condensation of charged p mesons 

The basic idea of our work is as follows. Consider 
a charged relativistic spin-s particle moving in a back- 
ground of an external magnetic field. Without loss of 
generality we assume that the magnetic field B cxt = 
(0,0, Bext) is directed along the z-axis, B cyA > and 
we consider spatially uniform and time-independent ex- 
ternal fields only. The energy levels e of the free particle 
of the mass m in the magnetic field are characterized by 
three parameters: the nonnegativc integer n ^ 0, the 
spin projection on the field's axis s z — — s, . . . , s, and the 
particle momentum along the field's axis, p z : 



< a >,)=pj + (2n-2«, + l)eB a 



(1) 



In this work we consider the charged particles, pions 
(s = 0) and the vector particles, p-mesons (s = 1), for 
reasons that will be clear later. For a moment, we as- 
sume that these particles are free, so that their (squared) 
minimal effective masses, corresponding to lowest energy 
states with p z = 0, are, respectively: 



m v ±(B ext ) 
m 2 p± (B cxt ) 



± + eB ext 
eB 



'ext 



(2) 
(3) 



The zero-field vacuum masses of the and p ± mesons 
are, respectively [14], 



139.6 MeV, 



775.5 MeV 



(4) 



Equation ([3| implies that the lowest energy of the 
charged p-meson in the external magnetic field may be- 
come purely imaginary if the magnetic field exceeds the 
following critical value 



Br = 



10 lb Tesla , 



(5) 



This observation indicates that the strong magnetic field 
(I?ext > B c ) makes the QCD vacuum unstable towards 
condensation of the charged p mesons. This new QCD ef- 
fect is very similar to the magnetic-field-induced conden- 
sation of the VT-bosons which was predicted by Ambj0rn 
and Olesen [12l 03] . The behavior of the lowest mass ^ 
of the charged p* meson in the region ^ B cxt ^ B c is 
shown in Fig. [ljby the solid line. 
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FIG. 1. Masses of the lowest p-meson eigenstates and of the 
products of their dominant decay modes as functions of the 
external magnetic field B = Boxt- The left red point and 
the middle blue point mark the onsets of the "7r ± -stability 
regions" for the neutral (111 and charged H p mesons, re- 



fer the neutral (111 and charged 
spectively. The right green point marks the critical field B c 
which corresponds to the onset of the condensation |5|. 

The subtle point of Eq. g [and of Eq. for s = 1 
as well] is that the gyromagnetic ratio of the vector p^ 1 
meson is set to be g = 2. In fact, this g- factor is "anoma- 
lously" large compared to the standard gyromagnetic ra- 
tio <7 m i n = 1 of a charged vector particle which is min- 
imally coupled to electromagnetism. Notice, that it is 
the anomalous gyromagnetic ratio g\\r = 2 which drives 
the condensation of the W bosons in the strong mag- 
netic field [T2l fTB"] . The large g- factor for W boson is 
a direct consequence of the non-Abelian nature of elec- 
troweak gauge group. 

As for the p mesons, the electrodynamics of these par- 
ticles has also elements of a non-Abelian structure which 
is visible in phenomenological Lagrangians [T5HT7] . The 
vector dominance hypothesis [T5] as well as the QCD sum 
rules [12] point out that the g- factor of the p- mesons is 2. 
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We discuss the quantum electrodynamics for these vector 
mesons in more details in Sec. Mil 

B. Larger lifetime of charged and neutral p mesons 

In the absence of the external magnetic field both the 
charged and neutral p mesons are very unstable particles 
characterized by the mean lifetime t p ps 4.5 x 10~ 24 s rj 
1.35 fm/c which corresponds to the full width [T4"] 

r p ^ an = 149.1 ±0.8 MeV. (6) 

Thus, one may incorrectly conclude that if even the p- 
meson condensate is formed at the strong magnetic fields, 
then it will anyway be unstable due to the intrinsic insta- 
bility of the p mesons themselves. Below we show that 
this statement is incorrect. 



1. Charged vector mesons 

Consider first the charged vector mesons. All known 
decays of the p£ mesons are going via the modes [14] : 

P ± ^tt ± X, X = tt°, Tj, 7, iririr. (7) 

The fraction of the primary decay mode, X — ir°, is 
greater than 99%. 

As the strength of the background magnetic field in- 
creases, the product of the decay, the charged pion [which 
is always created in the known decay modes of the p 
mesons Q] becomes heavier ^ while the decaying parti- 
cle, the lowest state of the p ± meson, becomes lighter 
Obviously, at a certain magnetic field B„± the masses 
of the initial and final states in the dominant channel, 
pn — > ^ir , should become equal, 

m P ± ( B P ± ) = m 7r± (B p ± ) + m w o , (8) 

and the fast strong decays ^ of the charged p mesons 
should eventually become impossible due to obvious kine- 
matical reasons. The strength of this "^-stabilizing" 
field is approximately three times weakeiQ compared to 
the critical field of the p condensation ([5]), 

B p ± = ~[ml-ml-m v (ml + 2m 2 p )^] ~ 0.36B C . (9) 

The left and right hand sides of Eq. ([8]) are shown by the 
solid and dot-dashed lines in Fig. [lp The point of the 
intersection of these lines gives us the critical field ([9]). 

At B > B p ± the charged p mesons may in principle de- 
cay via other slower (and undetected so far) channels that 
avoid fast gluon-mediated ir^ production. On the other 



Here and below we always neglect the difference between the 
masses of the charged tt^ 1 and p^ mesons, and their neutral 
counterparts, it and p°, respectively. 



side, the QCD string (which confines the quarks and an- 
tiquarks into mesons and baryons) is partially stabilized 
by the external magnetic field [2D]. Thus, in the suffi- 
ciently strong magnetic field the allowed modes of the 
decays of the charged p mesons should be much slower. 
As the result, the lifetime of the p ± -mesons should be 
much longer compared to the lifetime of these particles 
in the absence of the external magnetic field. 

One can also make a qualitative prediction for the be- 
havior of spectral function of the charged p meson in the 
strong magnetic field. Expected behavior of the lowest- 
mass peak is plotted in Fig.[2]as a function of an invariant 
mass. At zero magnetic field the p ± meson is seen as a 
broad resonance (the right peak in Fig. [2]) . As we switch 
on the background magnetic field, the single peak should 
split into multiple peaks corresponding to different levels 
of the charged vector particle (s = 1) in the external mag- 
netic field . The increase of the strength of the back- 
ground magnetic field leads to the kinematical suppres- 
sion of the p-meson decay modes, and, consequently, to a 
narrower lowest-mass peak in the corresponding spectral 
function (the peak in the middle of Fig. [2j. At B ^ B c , 
the onset of the condensation of the p-mesons occurs. 
This effect can be seen as appearance of a singularity of 
the (5-function-type located at the zero invariant mass 
(the left peak in Fig.[2|. 




Inv. Mass [GeV] 

FIG. 2. Prediction: a qualitative effect of the external mag- 
netic field on the lowest-mass peak in the spectral function of 
the p ± mesons vs. the invariant mass. As the magnetic field 
strength B = B oxt increases, the broad peak in the unstable 
low-B phase (right) turns into a much narrower peak in the 
7r ± -stable phase (middle). At the onset of the condensation 
of the p-mesons the peak transforms into the <5-function-like 
singularity located at the vanishing invariant mass. Features 
at higher invariant masses are not shown. 



2. Neutral vector mesons 

Similarly to its charged counterpart, the neutral p° me- 
son should also be ir ± -stabilized in a sufficiently strong 
magnetic field background. The primary channel of the 
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p° decay, p° —> tt + tt~ (it corresponds to more than 99% 
of the decays) becomes inoperative due to the same kine- 
matical reasons provided B cxt ^ B p o, where 



nipo(Bpo) = 2m 7r ±(B p o) . 



(10) 



The p° mass is expected to be practically independent of 
the magnetic fielcrl m p o(B) ~ m p o(B = 0), so that 



mi — 4:mz 

B a o = ~ 0.225 c 

p 4e 



(11) 



The left and right hand sides of Eq. (To] ) are shown by 
the dashed and dotted lines in Fig. [I] The intersection 
of these lines occurs at the critical field p"T] ). 

In the absence of the external magnetic field the neu- 
tral p meson has also other decay channels which do not 
involve the production of the charged ir^ pions. Such 
decay modes, however, are much slower compared to the 
primary decay 7r ± -modes like p Q — > ir + ir~. For example, 
the most effective 7r ± dess decay of p° is p° — > 7r°7, with 
the width 



0.089 ± 0.012 MeV. 



(12) 



which is more than three orders of magnitude narrower 
compared to the full width Q. In this paper we are 
not discussing how the 7r ± -less decays are affected by 
the strong magnetic field. However, it is clear that 
the electromagnetically-driven decay channels should be 
slower compared to the strongly-mediated ones. Thus, 
there are good kinematical reasons to believe that the 
prolongation of the p-meson life - induced by the strong 
magnetic field background - should be substantial. 

As for the evolution of the p° peak in the spectral func- 
tion, we expect that the background magnetic field makes 
it narrower, while its position is largely unaffected by the 
external field. As it will be clear from the results reported 
below, at B cxt > B c we may expect an appearance of a 
singular peak at zero p°-meson mass due to (quite weak, 
though) condensation of the neutral p° mesons. 



3. Reversed decays and effects of chiral condensates 

The estimations of the values of the critical fields ^ , 
([9]), and (11) are obviously approximate, as one may ex- 



pect systematic corrections coming from other effects of 
the strong magnetic field on the mass spectrum of the 
mesons. For example, in our qualitative considerations 
we don't take into account effects of mixing of the p° me- 
son with the neutral u> and (p mesons. We also neglect 
influence of the magnetic field on the p mesons and pions 



2 Here we ignore a weak coupling of the magnetic field to the mag- 
netic dipole moment of the p° meson. This coupling makes the 
critical field slightly stronger. 



at the quark level. However, the latter effect may be es- 
timated, at least partially. Indeed, the background mag- 
netic field enhances the chiral symmetry breaking |21) . 
According to a leading order of the chiral perturbation 
theory [22] (confirmed by the results of the recent lattice 
simulations [23]) the chiral condensate S is a linearly in- 
creasing function of the strength of the external magnetic 
field B ext : 



](0)(l 



In 2 

327^/1 



eB n 



(13) 



where / ff = 92.4 MeV is the pion decay constant. At 
the critical fields ( [IT] ), (|9| and Q the corrections (13) to 
the chiral condensate are 3%, 6%, and 16% respectively. 
We expect that uncertainties in our estimations of the 
critical values ( 11 ), Q and ^ may be of the same scale 
at least. 

One should also note that our considerations imply 
that at B > B p ± (B > B p o) the charged pions may de- 
cay into the charged (neutral) p-mesons. The statement, 
that the presence of the strong enough magnetic field 
interchanges the decaying and created particles, should 
not be disappointing. For example, it is known that the 
magnetic field may reverse the /3-decay of the neutron be- 
cause at the background magnetic fields with the strength 
greater than 5 • 10 14 T w O.lm^/e the proton becomes 
heavier than the neutron. As a consequence of this ef- 
fect, the proton may decay into the neutron by positron 
emission [2"4] . 

Summarizing this section, the charged and neutral p 
mesons are very unstable particles provided the magnetic 
field is weaker than the critical values B p ± , Eq. (|9| , and 
B p o, Eq. (11), respectively (Fig.[lJ. We expect, however, 
that as the external field becomes stronger than these 
critical values, the corresponding p mesons get stabilized 
with respect to the vast majority of the strong decays 
which are going via the production of the mesons (we 
call these regions of the magnetic field intensities as the 
"tt* -stable" phases both for the charged and neutral p 
mesons). If the background field surpasses the critical 
value (|5]), B cxt > B c , the condensation of the charged p^- 
mesons should occur. Below we show that at the same 
point -B xt = B c the neutral p° mesons may simultane- 
ously form an inhomogeneous superfluid. 



III. ELECTRODYNAMICS OF p MESONS 

A. The DSGS Lagrangian 

The self-consistent quantum electrodynamics for the 
p mesons was recently constructed by Djukanovic, 
Schindler, Gegelia and Scherer (DSGS) in Ref. [T^] start- 
ing from an effective Lagrangian for vector mesons devel- 
oped by Weinberg [15] long ago. The chiral, Lorentz and 
discrete symmetries of the Weinberg Lagrangian were ex- 
tended to the Maxwellian U(l) sector by adding all al- 
lowed interactions with electromagnetic fields. In terms 
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of the renormalized fields the bosonic part of the DSGS 
Lagrangian reads as follows [16] : 



r — — F F^ v — 



1 



(14) 



25s 



where A p is the photon field, = (pft — ip^)/y/2 

and = p^jp are, respectively, the fields of the (nega- 
tively) charged and neutral vector mesons^] characterized 
by the mass m p . The DSGS Lagrangian possesses the 
U(l) gauge invariance: 



pl°\x) -> p$>(x), 
p M (x) -> e-W P(I (x), 
A M (x) -> A A1 (x) + 9 AI w(a;). 



f/(l) 



(15) 



The tensor quantities in ( 14 1 are 



Fn V — dnA v — d v A, 



f(0) _ g (o) _ q JO) 



P^u = f$ - i9s(plpu - Pupl) 
PiLv = D^p v - D u p ti , 

and the covariant derivative is 



(16a) 
(16b) 

(16c) 
(16d) 

(17) 

= Pt 

meson fields carry the electric charges — e and +e respec- 
tively (here e = \e\ is the elementary electric charge). 
The coupling constant g s can be estimated [TBI H3 from 
the Kawarabayashi-Suzuki-Riadzuddin-Fayyazuddin re- 
lation [251: 



D„ 



i9spf 



ieA,, 



Equation ( 16 ) indicates that /?„ = p M and /v 



9s 



9p7T7T 



= 5. 



(18) 



so that g s >• e = v^Tr^cm. w 0.303. 

The most important fact for us is that the last term 



of the DSGS Lagrangian (14) describes a nonminimal 



coupling of the p mesons to the electromagnetic field. 
This term has two parts, 



5C = 5C {a) + 5C 



ch 



(19) 



SC^ = ^-(d,p^-d^)F^, (20) 

-ch _ „•„„ 3 f p^v _ 



5£ c 



(21) 



where the first part SC^ corresponds to the coupling of 
the electromagnetic field to the magnetic dipole moment 
of the p° meson, while the second part S£ ch describes 
the nonminimal coupling of the charged p -mesons to 



the electromagnetic field. The presence of the former 
may lead to an instability of the vacuum of the neutral 
vector particles (p^ mesons in our case) [55], while the 
latter implies the anomalous gyromagnetic ratio (g = 2) 
of the charged p ± mesons, so that the magnetic dipole 
moment of the p^ mesons is 



P P ± = ±2 • s. 

2m p 



(22) 



(here s is the meson's spin). It is the coupling (21) that 
plays a dominant effect in our paper while the interac- 
tion ( 20 ) makes a sublcading contribution. 



As we have already discussed in Section [II A[ spin-one 
particles with the gyromagnetic ratio g = 2 in strong 
enough external magnetic field should experience a tachy- 
onic instability towards development of a Bose-Einstein 
condensate. Since the condensed particles are charged, 
the condensate should be superconducting, and this fact 
is our central observation which is discussed in details 
below. 



B. Equations of motion 



A variation of the DSGS Lagrangian (14) with respect 
to the electromagnetic potential A p provides us with the 
Maxwell-type equation of motion, 

d v F vtl = - J„ , (23) 

where the electric current J M contains two contributions, 

(24) 

coming from the charged and neutral mesons, 

Jf = ie[p^p Vil - pfp^ + d\pl Ptl - p\p v )\ (25a) 
= ie[(D^ Pl/ -p^D^ 

+d»Up, - plp v ) + plD v Pil - (D v p^ Pv \ , 



9s 



(25b) 



respectively. The currents (25) are separately conserved: 
d» J M = d" Jf = = . (26) 



A variation of the DSGS Lagrangian (14) with respect 



to the field p^ gives us the second equation of motion, 

Pp® + mlpf--d"F Vil -ig s {p\ v p v - p^)=0 . (27) 
9s 



It can be rewritten as follows [we used (23), (24), (25)]: 



(Pd v + m% 0) )pM - ^,,(0) _ il Jf = o , (28) 



3 We denote the field of the neutral meson as (x) in order to 
discriminate it from the timelike component p° (x) of the charged 



so that Eq. ( 26 ) gives us 



(29) 
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Equation ( 28 ) provides us with the mass of the neutral 



p(o) m eson: 



too = m pi0) = to p 



(30) 



Using Eqs. (|24h , (25 1 and (28) one can get a well-known 



relation (emerged originally in the scope of vector dom- 
inance models long time ago [TH]) between the electro- 
magnetic current J p and the neutral meson field 



(0). 



Ju 



M 0) 



(31) 



(notice that in our notations e = |e| > 0). 
The third equation of motion is 

D"p vlt + m 2 pPp + i(g s p$ - ei^K = . (32) 
Using the identity [D p , D v ] = i(g s fji J - eF pv ), one gets 
(D a D a +m 2 p )g pu -D p D u 
+i{9sP$ + 9sfj?J - 2eF flv )] p v = . (33) 



Equations (|32|) and (|3l|) imply that 

D„ 



iff 



em 



■2 



f = . (34) 



The linear part of Eq. ( 33 1 gives us the mass of the 
charged p^ meson, 

m p ± = m p . (35) 

The neutral vector p^ meson is heavier comp ared t o its 
charged counterpart p ± . According to Eqs. (fig]), (30) 



and ( 35 1 , the difference in the masses is very small |16j . 

47TQ e . m j2 



5m p = toq — to„± 



1 MeV . (36) 



This mass difference is consistent with the available ex- 
perimental bounds 14J. 



IV. EXAMPLE: GINZBURG-LANDAU MODEL 

In the next Section [V] we analyze the condensation 
of the p mesons in the strong magnetic field, starting 
from the phenomenological field-theoretical DSGS La- 
grangian (14). However, before going into the details 



of the p condensation in QCD, it is very useful to discuss 
a few basic properties of conventional superconductivity 
in the condensed matter physics. Below we concentrate 
on the Ginzburg-Landau (GL) model which provides us 
with a simplest phenomenological description of the su- 
perconductivity. 



A. The Ginzburg-Landau Lagrangian 

The relativistic version of the GL Lagrangian for a su- 
perconductor is: 

£gl - -\f^F^ + (S) M $)*S)^ - A(|$| 2 - v 2 ) 2 ,(37) 

where T)^ = d p — ieA p is the covariant derivative and $ is 
the complex scalar field carrying a unit|^] electric charge e. 

The ground state of the model ( [37] ) is characterized by 
the homogeneous condensate of the scalar field, $o = (&) 
with |$q| = V- I n the condensed state the mass of the 
scalar excitation, <5$ = $ — $o, and the mass of the 
photon field A p are, respectively, as follows: 

to$ = v4\t) , rriA = V%er] . (38) 
The classical equations of motion of the GL model are 



S^i- + 2A(|$| 2 - 77 2 )$ = . 

d v F v » + J£ L = . 



where the electric current is 



7^ - -?> 

GL — (e 



(39) 
(40) 

(41) 



B. Destructive role of magnetic field 

The superconducting state in the GL model is com- 
pletely destroyed ($ = 0) in a background of the strong 
magnetic field B cxt , if the strength of the field exceeds 
the critical value 



B? h = 



~2e ~ 



2A 



-T 



(42) 



Let us assume that B cxt = -F12 is the only nonvanishing 
component of the field strength tensor. Consider the case 
when the uniform time-independent mag net ic field B cxt 
is slightly smaller than the critical value (42), B < B~ L , 
so that 



Bp-xf 
1 — <s- l 

Bc GL 

Then the condensate is very small 



(43) 



(44) 



and Eq. ( [39] ) can be linearized: 

{(©! - i2) 3 )(©i + *S 2 ) + e[B c - B(x)]}$ = , (45) 

where B(x) is the field inside the superconductor (here 
we consider static and z-independent solutions which cor- 
respond to a lowest energy of the system). In the vicinity 



4 Without loss of generality, it is convenient to consider the singly- 
charged bosons 4> instead of the usual doubly-charged bosons. 
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B c , so that Eq. (45) re- 



of the critical field B ~ B e - 
duces to the following equation for the condensate $: 



£»$ ~ with 23 = 33 1 + iD 2 ■ 



(46) 



The magnetic field destroys the superconductivity in 
the ordinary superconductor. On the contrary, we show 
below that a strong enough magnetic field should in- 
duce the superconductivity of the charged p mesons in 
the QCD vacuum. 



C. Abrikosov lattice of vortices in mixed state 



The GL model (37) admits a topological stringlike 



solution to the classical equations of motion (39) and 
(40), which is known as the Abrikosov vortex [27 . The 



Abrikosov vortices are formed when the superconductors 
are subjected to external magnetic fields. 

A single Abrikosov vortex carries the quantized mag- 
netic flux (remember that we consider the condensed 
bosons $ which carry the electric charge e and not 2e): 



J d 2 x±B(x±) 



2ir 

e 



(47) 



where the integral of the vortex magnetic field B is taken 
over the two-dimensional coordinates x± = (^1,^2) of 
the plane which is transverse to the infinitely-long, strait 
and static vortex. In the original solution, the scalar field 
of the unit-flux vortex is singular at the vortex center, 

(48) 



$(2j_) oc |xj_|e JV = x% + ix2 



where ip is the azimuthal angle in the transverse plane, 
and \x± 



is the distance from the vortex center. Equa- 
tion (48) corresponds to small \x±\: m$\x±\ <C 1 and 
mA\xT\-^ 1. 

In a type-II superconductor [in which m$ > rriA, or, 
according to (38), 2A > e 2 ] the Abrikosov vortices re- 
pel each other. If the external field is strong enough 
[but lower than the critical value (42)] then multiple 



Abrikosov vortices are created. Due to the mutual repul- 
sion, the vortices arrange themselves in a regular struc- 
ture known as the Abrikosov lattice [281 ES] . Since the 
normal (nonsupcrconducting) phase is restored inside the 
vortices, the Abrikosov lattice corresponds to a "mixed 
state" of the superconductor, in which both normal and 
superconducting states of matter are present. 

There are various types of the Abrikosov lattices which 
are characterized by different energies (25] ■ The stable 
lattice corresponds to a minimal energy of the system. If 
the magnetic field B cxt approaches the critical magnetic 
field (42) from below, then the simplest lattice type is 



given by the square lattice solution of Eq. ( 46 ) 



$(2:1,22) = C exp< 



(e-B ext ) 5 



(49) 



J 2 , o Xl + 1X2 
2 exp< — irn + Zirn — — 



In this equation the parameter Co is independent of the 
transverse coordinates x±. The inter- vortex distance Lb 
is expressed via the magnetic length £ B , 



2n£ f 



1 



VeB c 



(50) 



The area of the elementary square cell (i.e. of a cell which 
contains one Abrikosov vortex) is L B = 2tv£b- The abso- 
lute value of the condensate, |$(xj_)|, has a square sym- 
metry in the solution ( 49 ) and the vortices are located at 



the sites of the square lattice, 



Xj 

L B 



Hi € Z , 



1,2. 



(51) 



In this case the distance between the vortex centers is Lb ■ 



At the points (51) the condensate $(21,2:2) vanishes 
exactly and in the vicinity of these points the scalar 
field ((491) follows the behavior of Eq. (|48l. 



As we will see below, the pure superconducting state 
cannot be formed in the p meson superconductor con- 
trary to the ordinary superconductor. Instead, the 
Abrikosov lattice state is created. 



D. Homogeneous isotropic superconductivity 

Let us now apply a very weak external electromagnetic 
field to the superconductor. Neglecting effect of the ex- 
ternal field on the condensate $0, one gets from (41 ): 



^GL 



(52) 



where is given in Eq. (38). Setting p=0 and v=i in 



Eq. (52) one gets the first London relation for a locally 



neutral [Jo (2) = 0] superconductor 
9Jgjj _ 2 p 



(53) 



where E % = — F° l is the time-independent and uniform 
electric field. Equation ( 53 ) implies a linear growth of the 



electric current in external electric field, thus indicating a 
vanishing electric resistance of the superconducting state. 

In the long-wavelength limit, |<f| — > 0, the weak electric 
field E(x, t) = EQe %lyX '^~^ induces the local current 



J k (x,t;uj) = ^2cr k i(uj)Ei(x,t) 



(54) 



k=l 



where o~ki — Re <Jki + i Im o~ki is the complex electric con- 
ductivity. The London equation ( 53 ) indicates that 



Cfei(w) 



'kl 



(55) 



where the first contribution is a singular isotropic part 
associated with the superconducting state: 



2i 

1TUJ 



>kl 



(56) 



The regular part cr rcg accounts for all other (nonsuper- 
conducting) contributions to the conductivity. 

It is clear that the superconductivity described by 
Eq. ( 53 ) is homogeneous (it is independent of the spa- 
tial coordinate) and isotropic (it is independent of the 
direction). On the contrary, we will see below that a 
strong enough magnetic field induces inhomogeneous and 
anisotropic superconductivity of the charged p mesons in 
the QCD vacuum. 



E. Meissner effect 



The spatial components of Eq. ( 52 ) give us the second 
London relation: 



d x i/gl — -m%B , 



(57) 



so that in the absence of the external electric field (E = 0) 
one of the Maxwell equations (40), Jql = d x B, implies 



(-A- 



JB = 0. 



(58) 



This equation indicates that the photon inside the super- 
conductor acquires the mass wla, Eq. (38). Consequently, 



the superconductor tends to expel the external magnetic 
field ("the Meissner effect"). Physically, the Meissner 
effect is realized because the external magnetic field in- 
duces the circulating superconducting currents (57) in- 



side the superconductor. These currents, in turn, screen 
the external magnetic field since they induce their own 
magnetic field which is opposite to the external one (here 
we always assume that B ext < B c ). 

A weak magnetic field which is parallel to the boundary 
of the superconductor is always screened inside the bulk 
of the superconductor. The perpendicular magnetic field 
may penetrate the superconductor and create a mixed 
phase of the Abrikosov vortices. 

As we will see below, the second London equation 



( 57 ) is not realized in the superconducting phase of the 



QCD vacuum contrary to the conventional superconduc- 
tor. Consequently, the Meissner effect cannot be formu- 
lated in a selfconsistent way in the suggested supercon- 
ducting phase of QCD. 



V. CONDENSATION OF p MESONS 
A. Homogeneous approximation 



The energy density of the DSGS model (14) is 
1 



:F 2 



u o) ) 2 +(prr (59) 



„( n h 2 ' 



l 



+ m l(plpo + pi pi) - —FoiPoi 
9s 



were is the energy-momentum tensor, 
rp _ 9 _^C_ _ r 



(60) 



In order to understand the phase structure of the p 
mesons in the background magnetic field, it is useful to 
study first the homogeneous field approximation. To this 
end we ignore the kinetic terms d^pffl = and D^p v = 



in Eq. ( 59 1 . The remaining (potential) part of the energy 
density in the external uniform magnetic field B ext is: 



eoip^pW) = ~-B. 



+ -fW^-pW)] 2 (61) 



Tfl 

J r"J-^ (p\p-2 - iPl) + Jip}?)' 



™ 2 p pIp» 



where the sums over silent indices are written in the 
Euclidean metric, O 2 = X}u=o^V^V- We remind that 
we always take B cxt = -F12 > 0, and in this Section 
Foi = F 3 i = 0. 

The ground state of the model can be found via the 
minimization of the potential energy (61 1 with respect to 



the meson fields. To this end we notice that the field of 
the neutral meson is vanishing at the ener gy m inimum, 
p^ — 0. Then, the quadratic part of Eq. (61) becomes 
as follows: 

4 2) (Pm) = ieB ext (p[p 2 - p\pi) + m 2 p p\ Pil 



a, 6=1 



piM ab p b + m 2 (plp + plp 3 ) . (62) 



The Lorentz components p\ and p 2 possess the non- 
diagonal mass matrix 



M = 



ieB. 



—ieBr. 



ext 
2 



(63) 



The eigenvalues p± and the corresponding eigenvectors 
p± of the mass matrix (63) are, respectively, as follows: 



p\ = m 2 ± e-Bext 



1 



P± 



V2 



(pi + ip 2 )- (64) 



The mass terms for p and p 3 components are diagonal in 
( 62 ) and their prefactors m 2 are unaltered by the external 
magnetic field. 



It is clear from Eq. ( 62 ) that in terms of the "longi- 



tudinal" components po and P3, the ground state of the 
model corresponds to pg = p^ = at any value of the 
magnetic field. We express the transverse components 



P1.2 via the eigenvalues and eigenvectors (64) of the mass 



matrix (63), and then we get for (the potential part of) 



the energy density ( 59 ) the following expression: 



£ o(p + ,P-) = ^ c 2 xt + f (\ P+ \ 2 -\p-\ 2 f 

+/4| P+ | 2 + M 2 _|p_| 2 . (65) 
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Since p\ > regardless of the value of the magnetic field 
B cxt , the ground state corresponds to p + = 0. In turn, 
this means that p_ = y/2p and 



should exhibit spatially anisotropic properties due to spa- 



Pi 



"«P2 



PO = P3 = , 



(66) 



where p is a scalar complex field. In terms of the new 
field p the energy density ( 65 1 takes the simple form: 



2 

ext. 



tially anisotropic condensate ( 66 ) . This issue will be dis- 
cussed in details later. 

Note that in the presence of the background mag- 
netic field B cxt the rotational group SO(3) ro t is explic- 
itly broken to its O(2) rot subgroup generated by rota- 
tions around the axis of the magnetic field. In the homo- 
geneous approximation, the ground state ( 66 ) is trans- 



+ 2 (™p - e-Bcxt) |p| 2 + 2g 2 s |p| 4 ■ (67) formed under the global O(2f rot rotations aVfollows 



Thus, we get the familiar Mexican-hat potential which 
describes various spontaneously broken systems. In par- 
ticular, the same potential appears in the GL model of 
superconductivity ( 37 ) . 



The ground state of the model (67) depends on the 



value of the external magnetic field: if the field strength 
is weaker than the critical value B c = m 2 p /e, Eq. (pi, 
then the potential is trivial, while if -B cx t > B c then we 
get a nontrivial ground state: 



IpIo 



'e(B axt -B c 



B C xt ^ B c 
Bo^t < B c 



(68) 



(the subscript "0" in \p\ indicates that we consider the 
homogeneous-field approximation). In Fig. [3] we plot the 
behavior of the condensate ( 68 1 as the function of the 
external magnetic field i? GXt . The value of the condensate 
follows a typical behavior of an order parameter for a 
second order phase transition at B cxt = B c . 




eB, GeV 2 

FIG. 3. The condensate \pas\ of the charged p mesons as a 
function of the external magnetic field B = B oxt at the ground 
state. This single curve describes both the uniform conden- 
sate | Pas | = |p|o in the homogeneous approximation ( |68| l and 
the mean-cell value \pas\ = IpasI^ of the inhomogeneous 
condensate ( |96[ ) in the weak-amplitude approximation. 

In Fig. [3] the subscript "AS" in \pas\ stands for the 
"anisotropic superconductor" . Indeed, the scalar field 
p(x) enjoys the gauge symmetry ( 15 ) of its vector prede- 
cessor p^(x), 



U(l) e 



p{x) 



iu)(x) 



p(x) 



(69) 



0(2) rc 



p(x) -+ e^p(x) . 



(70) 



where ip is the azimuthal angle of the rotation in the 
transverse plane. Thus, the ground state (66 1 is invariant 



under a combination of the global transformation from 
the gauge group ( 69 1 and the global rotation around the 
field axis ( 70 ) provided the parameters of these transfor- 



mations are related ("locked") to each other as follows: 
u)(x) — ~tp. In analogy with the color superconductiv- 
ity [§] one can say that the ground state "locks" the resid- 
ual rotational symmetry with the electromagnetic gauge 
symmetry: 



U(l) c . m . x O(2) rot -► 1/(1) 

locked ■ 



(71) 



Below we will see that the inhomogeneities of the con- 
densate break the locked group ( 71 ) further to the group 



of discrete rotations of the vortex lattice. 



In the ground state ( 68 1 the potential energy ( 62 ) has 
the form: 



£o(M = Ho) 
eo(\ P \ = 0) 



«»-*) 



-Dext B c 

B C xt < B c 



(72) 



Obviously, for a strong magnetic field B cxt ^ B c , the 
condensed state has lower energy compared to the energy 
e (\p\ = 0) = Sg Xt /2 of the normal (noncondensed) state. 

Thus, we observed that the condensation of the p^- 
mesons in the QCD vacuum should be very different from 
the condensation of the Cooper pairs <& in the standard 
superconductor which is described by the phenomenolog- 
ical GL model (37). Indeed, in Section IV B we have il- 



lustrated the destructive role of the strong magnetic field 
on the conventional superconductivity. On the contrary, 
in this Section we have found that the strong enough 
magnetic field enforces the p-meson superconductivity. 



B. Two-dimensional equations of motion 

In order to study the properties of the emerged su- 
perconductor in more details we should definitely go be- 
yond the homogeneous approximation. The inhomoge- 
neous state can be treated with the full system of the 
3+1 dimensional equations of motion for the p-meson 

We notice, 



The formation of the nontrivial ground state p in the fields which were discussed in Section III B 
strong external magnetic field B cxt ^ B c breaks spon- 
taneously the gauge symmetry ( 15 ) and forms, conse- 
quently, a superconducting state. The superconductor 



however, that a wavefunction of the lowest energy state 
of a free particle in a uniform static magnetic field is 
independent on the coordinate X3 which is longitudinal 
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to the magnetic field. The dependence on the time co- 
ordinate xq comes as a trivial phase factor only. The 
Abrikosov lattice solution in the type-II superconductors 
is also known to be independent of xo and X3 coordinates 
(Section |IV C ). These well-known properties suggest us 
to concentrate on xq- and a^-independent solutions to 
the classical equation of motions for the p mesons. To 
this end we choose the complex coordinate z = x\ + 1x2 
where xj_ = (xi,X2) are the coordinates in the spatial 
plane which is transverse of the magnetic field axis. We 
define the complex variables 



O = O x + i0 2 , = Ox- i0 2 



(73) 



for the fields O = J (0) , J, p (0) , A, and for the derivative 
O = d. It is also convenient to introduce two covariant 
derivatives: 



D = Dx + iD 2 = £ + ig s p 



(0) 



£> = d - ieA. (74) 



For the sake of convenience we use below both x± and 
z notations interchangeably, so that the two-dimensional 
Laplacian, for example, can be written in the three dif- 
ferent ways: dd = dj_ = df 



dl 



Our homogeneous field analysis (Section V A I suggests 



that the charged currents should be chosen in the form 



Po = pa = , pi = -ip2 = p{z) 



(75) 



where p is a complex fielcpj 

The magnetic field ( 16a ) and the field-strength of the 
neutral vector bosons ( |16b ) are as follows 

(76) 
(77) 



F 12 = lm(BA) = B{z) , 



Notice, that despite the external magnetic field B cxt is 
assumed to be uniform, the magnetic field (76) of the 



classical solution may be (and, in fact, will be) inho- 
mogeneous. The tensor quantities ( 16c ) and ( 16d ) take, 



respectively, the following form (we omit the argument z 
hereafter): 



pf 2 ] =C + 2g s \p\' 



P12 = iDp. 



(78) 



The charged and neutral components of the cur- 
rent (25) become simple expressions, respectively: 



ch 



2ie(p^Dp + d\p\ 2 ) 



=i~dC. (79) 
9s 



The conservation law for the charged current (26), 
Im.{d\p^Dp — } = 0, is satisfied automatically 

due to relation pIL 



Dp = 0, 



(80) 



5 In a strong field limit one can show that due to presence of inho- 
mogeneities the ansatz ( |75| may be generalized : pi = p(z)+£ (z), 
P2 = i[p{z) — In our analysis we ignore the subleading field 

£ because its amplitude is suppressed by the factor e/g s <C 1. 



[we also used the identity 9|,o| 2 = p^ Dp + (Dp)^ p]. 
Equations (23), (28), (32) reduce, respectively, to 



g s dB + iemlp w = 0, 
2 5s 2 |p| 2 )p(°)-2 i5s 9H 2 =0, 
?)]P = 0. 



-dd 



-DD + 2(g s C-eB + 2g 2 s \p\ 2 



(81) 
(82) 
(83) 



Equation ([77]) along with the conservation law (|29j), 
He(dp(°)) — 0, lead to a simple expression for the trans- 
verse component of the field tensor ( 16b ) of the neutral 
mesons: 



C = -idpM 



(84) 



C. Inhomogeneous condensate of small amplitude 

1. Linearized equations of motion 



The classical equations of motions ( 80 1— ( 84 ) comprise a 



complicated system of equations which is difficult to solve 
analytically due to the nonlinearities. However, following 
our discussion for the GL model (Section IV B), let us 



assume that the amplitude of the condensate p is very 
small. Then, the equations of motion can be linearized 
and a leading analytical solution can be obtained. The 
condensate p should be small if the background magnetic 
field B ext exceeds slightly the critical value B c , Eq. ([5]). 
Concretely, for i? GX t *J B c we consider the condition, 



2 3 2 |p| 2 «m 2 



B r 



' 



or 



- K 1 . 



(85) 



These relations are analogous to, respectively, weak- 



condensate conditions (44) and (43) in the GL model 



of superconductivity. We show below that the first and 



the second relations in (85) are, in fact, equivalent 



Notice, that Eq. (80) coincides with Eq. (46) for the or- 



der parameter for the ordinary superconductivity in the 



GL model (37) provided that the external magnetic field 



is close to the critical field ( 42 ) of this model. Therefore 



we should expect emergence of an analogue of the vortex 
lattice ( 49 ) in the p system ( 14 ) similarly to the appear- 



ance of the Abrikosov lattice (49 1 in the GL model. Thus, 



the condensate of the p mesons in the external magnetic 
field should definitely be inhomogeneous. Following the 
classic example }29| , we consider below the simplest case 
of the square lattice with the elementary length ( 50 ) . 



In the weak-condensate regime ( 85 ) we can work in 



the leading order in terms of the condensate pas (higher- 
order corrections are always omitted below). Then the 
equation of motion ( 82 ) gives the following relation: 



where 



-d\ + ml 



2ig s 



-3\ 



!<9|PAS| 



2tt 



K (m\x±\ 



(86) 



(87) 
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is the two-dimensional Euclidean propagator of a scalar 
particle with the mass mo and Kq is a modified Bessel 
function (remember that the subscript "AS" stands for 
the "anisotropic superconductor" solution). 



and its solution becomes as follows 



B AS (xj_) = B c 



2emQ 



I PAS I 



2e(| Ms | 2 )^.(92) 



It is very important to notice that Eq. ( 86 ) relates the Here the last term 



condensate of the neutral p° mesons with the condensate 
of the charged p^- mesons. Thus, if the have an inhomo- 
geneous condensate of the charged p ± mesons, then we 
automatically get the inhomogeneous condensate ( 86 1 of 
the neutral p° mesons as well! This fact may indicate, 
that the superconductivity of the p mesons may induce 
the superfluidity of the p° mesons. Notice that a relation 
between the superfluidity of the neutral p^ mesons and 
the superconductivity of the charged mesons may be 
guessed from the fact of the vector dominance, Eq. (31 1. 



We interpret the nonzero condensate ( 86 ) as a "super- 
fluid" because of the complex nature of the field p(°\ 
Moreover, if, for a moment, we assume that this field 
is homogeneous (i.e., coordinate- independent) then rota- 
tions of the system around the magnetic field axis ( 70 ) 



would transform it as a usual complex field in simplest 
bosonic theories of superfluidity [55], p^ — > e lv p^ . The 
inhomogeneities of the condensate (86) break sponta- 



neously this global group down to a discrete group of 
the rotations of the vortex lattice. 

We would like also to notice an important role of the 
inhomogeneities in the charged p ± condensate for the su- 
perfluidity. In Section |"V A| we have seen that the homo- 
geneous condensate of the charged p ± mesons alone is 
unable to induce the superfluidity of the p' ' mesons: a 
uniform nonzero expectation value of p does not imply 
p^ =/= 0. However, the inhomogeneous charged conden- 
sate of p^- automatically induces the inhomogeneous neu- 
tral condensate of p^ as one can see from the presence 
of the derivative d in the numerator in the right hand 
side of Eq. ((86). 



The transverse component of the stren gth tensor ( 77 ) 
of the neutral p° mesons is given by Eq. ( 84 ) : 



C A s(x±) = -idpfl = 2 9s _q 2 ~ 



y ± -r ,;* 



Due to the identity, 



d 2 

d 2 x ± o9 ± 2 (x±-y±) = 0, (89) 



-d 2 



the total "flux" of the neutral p° mesons through the 
transverse plane is always zero, 



(\PA S \ 2 ) A = J Y d 2 y ± \p AS (y ± )\ 2 . 

B J A 



(93) 



is "the mean-cell value" of the condensate squared |/?as| 2 - 
Due to the identity, 



d 2 x± 



-d\ + ml 



(xj_ - yj_) = 1 . 



(94) 



the last term in Eq. ( 92 ) guarantees the conservation of 



the net magnetic flux through each elementary cell A: 



f d 2 x ± B AS (x ± ) = f d 2 x ± B ext = L%B 

J A J A 



' xt -- — , (95) 

e 



here we have used Eq. ( 50 )] . The quantization of the 



magnetic flux ( 95 1 is similar to the quantization of the 
flux of the Abrikosov vortex (47). 



Finally, Eqs. (83 1, (90 1, (92) and (95) give us 



\PAS\ 




B cx t 5 s B c , 



Bext < B c 



'(96) 



for the mean value ( 93 1 of the condensate. The mean-cell 
value of the condensate (96) is shown in Fig. [3) Notice, 



that the mean-cell value of the condensate ( 96 ) coincides 



with the value of the uniform condensate (|68|) obtained 
in the homogeneous-field approximation. 

Equation (96) has a few interesting properties. Firstly, 



this equation represents a typical behavior of an order 
parameter. Secondly, Eq. (96) suggests that the phase 



transition, which separates the superconducting and the 
nonsuperconducting phases at B oxt = B c , is of a second 
order (as it is seen clearly in Fig. [3| . And thirdly, Eq. ( 96 ) 



proves the equivalence between the first and the second 
conditions of the weak-condensate regime, Eq. (85). 



Concluding this section we would like to stress that 
here we have introduced the new topological object, the 
" p vortex" , which is the vortex made of the supercon- 
ducting p^- mesons and superfluid p° mesons. This unit 
vortex-cell car ries the nonzero q uant ized flux of the mag- 
netic field ( 95 1 and zero /9°-flux ( 90 ) . The lattice of such 



vortices is a ground state of the superconductivity of the 
QCD vacuum at strong magnetic field. We discuss this 
lattice state in details in the next section. 



A 



d 2 ^± (/i ( 2 0) )as(^) = / C AS (x ± ) = 0, (90) 



where the integral is taken over a unit cell A of the peri- 
odic structure of the "p vortices" . 



Next, Eq. (81) gets simplified, 



B 



2em\ 



I PAS I 



(91) 



2. Inhomogeneous condensate: p-vortex lattice 



In the regime ( |85[ ) the degree of the inhomogeneity of 
the magnetic field SB(xj_) = B A s(x±) — B cxt in the su- 
perconducting state is extremely small. Indeed, accord- 
ing to Eq. (l92l, 



\SB\ ~ 2e\p\ 2 < 



r-Bcxt <§C B c 



(97) 
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Thus, the inhomogeneity of the magnetic field SB is 
supp ress ed both by the small amplitude of the conden- 
sate (85) and by the very small factor e 2 /g 2 =8.8xl0~ 3 . 
From Eqs. ( 88 ) and ( 92 ) one also finds that the stress ten- 



sor of the neutral bosons ( 77 ) is small compared to the 



magnetic field (76), |C| <C (e/g s )-B xt- 

Therefore we can set below B(x) ~ B CKt with the very 
good accuracy. Then, 

eB 

£~£ cxt = d-eA oxt =d+ — z (98) 



so that the solution of Eq. ( 80 ) is 



Pas (z) 



= P 4 



H AS (z/L B ) 



(99) 



where Has(z) is arbitrary analytic function of the ar- 
gument z and the inter- vortex distance Lb is given in 
Eq. ( 50 1 . Following the known solution ( 49 ) in the con- 



ventional superconductivity [29], we choose the square 
form of the lattice cells. For such periodic structure one 
gets 



Ha&(z) 



< e(_B cxt — B c ) 



V2g 2 s 



+ 00 

E 



-irn -\-2irnz 



,(ioo) 



where the prefactor was determined with the help of the 



expressions (99|) and (|l00|) 
We already 



normalization relation ( 96 ) supplemented by the explicit 



enow that the homogeneous condensate 



locks the rotational and gauge degrees of freedom ( 71 1 



The inhomogeneities in the condensate break the locked 
subgroup ( 71 ) further down to a discrete subgroup of the 



lattice rotations G]^ cked : 



f/(l)e.m. X O(2) rot -> U(l) locked ^ G 



lat 

locked ' 



(101) 



The discrete group G}o*kcd depends on the lattice struc- 
ture formed by the vortices. 

Similarly to the mixed state of the ordinary type-II su- 
perconductivity, the p-vortex centers are located at the 
points (51), where the condensate pas vanishes. In the 



vicinity of the p- vortex centers the condensate ( 99 ) fol- 
lows the typical Abrikosov- vortex behavior (48). How- 
ever, there are many essential dissimilarities between the 
vortex systems in the GL model and in the system of the 
condensed p mesons. 

In Fig. [4] we visualize four elementary lattice cells of 
the p vortex lattice in the transverse plane. We take 
the external magnetic field with the strength eB cxt = 
(800 Me V) 2 > eB c , so that the system is already in the 
superconducting state. The strength of the field satis- 
fies the weak-condensate condition (85). The magnetic 



length and the elementary distance between the vor- 
tices in the square vortex lattice are, respective!}^ (50), 



Note that the magnetic length £b is of the order of the size of the 



p meson itself, r p 



0.25 fm. Thus, at these magnetic fields 



is = 0.25 fm and Lb = 0.63 fm. The mean value of the 



condensate (96) of the p mesons is \(f>As\ — 23MeV. 
In Fig. [4] we plot various quantities that characterize the 
vortex: the amplitudes of the superconducting and super- 
fluid condensates, the excess of the magnetic field with re- 
spect to the external magnetic field and the field strength 
of the neutral meson field C. One can clearly see that: 

• The superconducting conde nsate p of the charged 
vortices p ± , Eqs. (99) and ( 100| , vanishes at the 
centers of the vortices (51), Fig. |4|a). In the vor- 
tex core the amplitude of the condensate \p\ is a 
linear function of the distance from the vortex cen- 
ter. This feature is similar to the behavior of the 
condensate near a typical Abrikosov vortex with a 
unit vorticity (48). 



• The superfluid condensate p(°\ Eq. (86), has a 
toothlike structure, Fig. |4]Jb). It vanishes at the 
locations of all local extrema of the superconduct- 
ing condensate including the centers of the vortices. 
The amplitude of the superfluid condensate is max- 
imal at the points of steepest behavior of the su- 
perconducting condensate, Fig.Qa). 



The magnetic field strength B, Eq. (92), takes 



its minimal values at the centers of the vortices, 
Fig.|4jc). The maxima of B are located outside the 
vortex cores. This feature contradicts our intuition: 
in the ordinary superconductivity the strength of 
the magnetic field is maximal at the center of the 
Abrikosov vortex. In fact, the p^- condensate has 
its own magnetic dipole moment due to the large, 
g = 2, gyromagnetic ratio of the p vortex. This 
dipole moment contributes only to the magnetic 
field outside the vortex cores, where the conden- 
sate of the p ± condensate is large, Fig. |4^a) . The 
electric current J, Eq. (24), is visualized in Fig. [5] 



• The strength of the neutral meson field C of the 
superfluid [Eq. ( [88] )] takes its maxima at the loca- 
tions of the p vortices, Fig. |4|d). Thus, the p vor- 
tices share this important property of the ordinary 
superfluid vortices as well. 

Summarizing, the vortex core expels both supercon- 
ducting and superfluid condensates of the charged and 
neutral p mesons, respectively. The magnetic field takes 
its maxima outside the vortices, while the strength of 
the superfluid (electrically neutral) field is peaked at the 
vortex centers. 



the p mesons should mutually overlap similarly to the overlap- 
ping Cooper pairs in the conventional superconductivity. How- 
ever, regarding the success of the phenomenological GL model of 
the superconductivity we don't question the applicability of the 
phenomenological DSGS model ^ I4[ ) in the strong-field regime. 
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FIG. 4. Four elementary cells of the p-vortex lattice in the 
plane x± — {x\,Xi). The plane is perpendicular to the ex- 
ternal magnetic field with the strength eB cxt — (800 MeV) 2 . 
From top to bottom: (a) the amplitude of the superconduct- 
ing condensate p, Eqs. (99 1 and jlOOl; (b) the amplitude of 
the superfluid condensate p^°\ Eq. ( |86| ; (c) the excess of the 
magnetic field SB(x ± ) = B(x ± )-B c ^ t , Eq. {92|; (d) the field 
strength C of the superfluid condensate p^°\ Eq. (881. 
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FIG. 5. The transverse components Ji and J2 of the electric 
current J, Eq. (241, in the transverse plane x± = (xi,X2)- 



Four elementary cells of the p-vortex lattice at the external 
magnetic field with e_B ex t = (800 MeV) 2 are shown. 



3. Anisotropic superconductivity 

The basic property of a superconductor is the absence 
of the resistivity. This feature is reflected, in particular, 
in the first London equation ( 53 ) in the GL model. 



There is a simple way to derive analogues of the Lon- 
don equations for the condensed state of p mesons in the 
external magnetic field. First, we notice that Eqs. (28) 
and (31 ) imply: 



(d a d a 



d[„J,,] = m n duJ', 



ch 



u \p. J v ] 



(102) 



Then, we take p — and v = 3 in Eq. (102) and use 



Eq. (25a) to get expressions for the p — 0, 3 components 



of the charged currents: 



j: h = 2 l e[p*D a p-(D a py P ] 



0,3. (103) 



Following the logic of the derivation of the London 
equation ( 53 ) in the Ginzburg-Landau approach (Sec- 
tion 



IVDj), one gets from Eqs. fl02| ) and ( [103 1 

dJ 3 (x ,x±) 



dx 



Ae 2 h 2 AS (x ± )E 3 



(104) 



where xj_ — (xi,x 2 ). The inhomogeneous quantity 



h 2 - 
"as — 



I Pas I 



(105) 



plays the role of the |$o| 2 condensate or m\/e 2 in the 
conventional London relation 1531). 



Equation (104) implies that the p-meson condensate 



exhibits the superconductivity phenomenon along the di- 
rection of the external magnetic field B: the electric cur- 
rent growing linearly with time if a weak external electric 
field is applied. 
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Notice, that due to the periodicity of the inhomoge- 
neous condensed state the mean-cell values of the squares 



of the "effective condensate" ( 105 ) and of the real con- 



densate ( 93 ) coincide identical 



(hls) A ^(\PAs\ 2 ) A . 



(106) 



Averaging Eq. ( 104 ) over an elementary square cell in 



transverse directions and using Eq. ( 96 1 we get the cell- 
averaged value of the electric current (J 3 )_4: 



di 



(Js)a 



2e 3 

2" (Bextt 



(107) 



quenched QCD vacuum which was found in lattice sim- 
ulation in Ref. [8] for weaker magnetic fields. An expla- 
nation of the anisotropy could be as follows: in a back- 
ground of a uniform magnetic field the electric charges 
may move along the axis of the magnetic field while the 
motion in the transverse direction is limited to the spatial 
size £b of the low Landau orbits (50 1. In a sufficiently 



strong magnetic field, and in absence of scattering of the 
charge carriers, the net transverse motion of the charges 
is suppressed contrary to the motion in the longitudinal 
direction. 



where B cxt > B c and we assumed, as usual, that the ex- 
ternal electric field E3 is a space-time independent quan- 
tity. The longitudinal (i.e., directed along B cxt ) super- 
conductivity sets in as the external field B cxt exceeds the 
critical value B c , Eq. 

It is easy to prove that the superconductivity phe- 
nomenon has an anisotropic nature: in the transverse 
(i.e., perpendicular to B cxt ) directions the superconduc- 
tivity is absent. In order to prove this fact let us ap- 
ply a weak spacetime-independent electric field E ext = 
(_E extj i, E CK t.2, 0) perpendicularly to the strong magnetic 
field background B cxt = (0,0, B cxt ). This electric field 
should test a possible transverse superconductivity of the 
p^-meson condensate which could also be created by the 
strong magnetic field. 

In order to show that the B cxt -transverse electric field 
does not create an accelerating electric current, we no- 
tice that an appropriate Lorentz boost may transform 
this system of the nonparallel E cxt and B cxt fields into 
the frame where the electric field is zero, E' ext = 0. Ob- 
viously, in the new frame there are no linearly growing 
electric currents, so that in the initial frame such runaway 
currents are absent as well and 



dJj(x ,xj_) 

at 



= 0, 



d_ 
at 



(Ji) A = 0, i= 1,2. (108) 



This argument does not work for the parallel electric and 
magnetic fields which were used to prove the longitu- 
dinal superconductivity (104 1. Indeed, in this case the 



scalar product (E ext ■ B ext ) cx E^apF^F^ is a Lorentz- 
invariant quantity which is insensitive to boosts and ro- 
tations. Thus, if -Eext II -E'ext then there is no frame where 
the external electric field E' ext is zero. 



Equations (107) and (108) imply that the (cell- 
averaged) electric conductivity ( 54 ) contains an 



anisotropic complex-valued contribution ( 55 ) which has 
a singular part at oj 



0: 



ire 

IT 



B„ 



B r . 



5{lv) 



2i 

7TCU 



SmSu, (109) 



where the index i — 3 corresponds to the direction of the 
external magnetic field B oxt . 

The anisotropy of the superconductivity is quite sim- 
ilar to the anisotropy of the "usual" conductivity of the 



4- Absence of longitudinal Meissner effect 

We have a very unusual situation: in our paper we sug- 
gest that in the QCD vacuum the strong magnetic field 
induces the superconductivity of p mesons, while all our 
experience in the condensed matter systems tells us that 
we should expect the opposite phenomenon ,28, 29] : the 
external magnetic field should destroy the superconduc- 
tivity due to the Meissner effect (Section IV E). In order 



to find a reason for this would-be inconsistency between 
the usual superconductor and the p-meson system let us 



apply the considerations of Section IV E to the p mesons. 

According to the Maxwell equations the electric cur- 
rents that could screen the external magnetic field -B ex t = 
(0,0, -B ex t) should circulate in the transverse x^-plane. 
In turns, the superconducting current in the transverse 
plane, Jas = -Tas 1 + ^as,2, can be related to the neu- 



tion (31): 



tral meson current ( 86 ) via the vector dominance rela- 



2iem\ 
-d\ + rr. 



,-5|pas| 2 -(110) 



Then in the system of the condensed p mesons, the ana- 
logue of the second London equation (57) for the longi- 



tudinal magnetic field can be written as follows (here we 
use the relation dd 



dl) 



(d x J AS ) 3 = Im(<9J A s) = 2em\ 



91 



-dl 



HpasI 2 (ih) 



The right hand side of this equation depends on the exter- 
nal magnetic field -B ex t via the superconducting density 

PAS, Eq. m. 



Equations (111), (99) and (100) provide us with an 



implicit expression for the curl of the screening currents. 
However, even without knowledge of the explicit form of 
these solutions one can show that these transverse cur- 
rents both screen and enhance the external magnetic field 
in such a way that the net effect in one elementary vor- 
tex cell is precisely zero. Indeed, let us integrate left and 
right hand sides of Eq. ( 111 ) over an elementary unit cell, 



take into account the periodicity of the solution ( 99 ) and 
use the following property 



d 2 x ± 



-dl 



r(a;_L - y±)=0. 



(112) 
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Thus, the cell-averaged right hand side of the second Lon- 
don equation (111 I for p mesons is zero, 



d 2 x_L (d x Jas)3 = , [condensed p^ mesons] .(113) 



while in the GL model the same procedure would give us 
the constant quantity in the right hand side of (|57j): 



d 2 Xl _ {d x j GL ), 



[GL model] . (114) 



field B" xt is to rotate the primary field B' ext . After sim- 
ple rotation of our coordinate around its origin we get 
a new field i? ex t = B' ext + B" xt so that the role of the 
additional "test" field is to rotate the directions of the 
p-vortices in the condensed state. Thus, the question of 
the (non) existence of the transverse Meissner effect can- 
not be formulated in a selfconsistent way. 



VI. CONCLUSIONS 



This fact simply means that in the state of the con- 
densed p mesons, the external magnetic field induce the 
transverse superconducting currents which are circulat- 
ing both clockwise and counterclockwise contrary. Con- 
sequently, the external magnetic field is enhanced in some 
regions of the transverse plane and it is suppressed in the 
other regions. Contrary to the ordinary superconductor, 
the net current circulation of the superconducting p cur- 
rents per a unit lattice cell is exactly zero (113), while 



in the ordinary superconductor the net circulation is a 
linearly growing function of the external magnetic field. 

Thus, we have found that the external magnetic field of 
any strength B cxt > B c does not experience the screen- 
ing inside the p-superconductor: the magnetic flux prop- 
agates freely inside the superconductor. The same state- 
ment is not true for the ordinary superconductor in the 
purely superconducting state: the magnetic field tries to 
avoid the superconductor (the Meissner effect). Thus, in 
a loose sense one can interpret the absence of the net 
circulating currents ( 113[ ) as the absence of the "longitu- 
dinal" Meissner effect. 

On the other hand, our system is very similar to the 
ordinary Abrikosov lattice in the mixed state of the type- 



II superconductor, Section IV C in the mixed state the 
magnetic field forms an inhomogeneous state and prop- 
agate though the superconductor, basically, in the cores 
of the Abrikosov vortices. In this case, however, the ex- 
ternal magnetic field must be bounded both from above 
and from below, contrary to our p-superconductivity in 
the QCD vacuum. 

One can try to address the question about the existence 
of the Meissner effect in the p-superconductor in a dif- 
ferent way. In the ordinary superconductivity the Meiss- 
ner effect is usually formulated as follows: if we apply a 
weak "test" magnetic field, say B' ext = (B' ext ,0, 0), along 
the boundary of a superconductor then this field will be 



screened inside the superconductor according to Eq. ( 58 1 , 
i.e. B(x 3 ) = {e-" lAX3 B' cxt ,0,0). This experiment, how- 
ever, is senseless in the case of the p-condensation be- 
cause this condensation is induced in the rotationally- 
invariant vacuum by the magnetic field itself. Indeed, 
assume that we have a combination of the two external 
magnetic fields: the strong field B' cxt , which induces the 
conductivity, and the additional weak field B" xt which is 
superimposed onto B' ext transversely {B' ext ■ B" xt ) = 0, in 
order to check the Meissner effect. Due to vacuum en- 
vironment it is clear that the sole role of the additional 



We argue that in a sufficiently strong background mag- 
netic field the QCD vacuum may undergo a spontaneous 
transition to a superconducting state via condensation 
of the charged p mesons. The critical strength of the 
magnetic field is given in Eq. ([5j . The superconductivity 
is understood in the usual electromagnetic sense. More- 
over, unlike the color superconductivity, the supercon- 
ducting QCD state is suggested to be formed in the cold 
vacuum, i.e. at zero temperature and at zero chemical 
potentials. Our vision of the phase diagram of the cold 
QCD vacuum in terms of the p-meson degrees of freedom 
is illustrated is Fig. [6] 
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FIG. 6. The expected phase diagram: the impact of strong 
external magnetic field B = B cxt on the p-meson degrees of 
freedom in the QCD vacuum. 

We have found the following basic properties of the 
superconducting state: 

1. The superconducting effect occurs because of the 
nonminimal coupling of the charged p mesons to 
the electromagnetic field. The strong magnetic field 
enhances the electromagnetic superconductivity of 
the QCD vacuum instead of destroying it. 

2. Due to simple kinematical reasons a strong enough 
magnetic field makes the lifetime of the p mesons 
much longer by closing the dominant decay chan- 
nels (p —> ^ir and p° — > ir + ir~) of the p mesons 
into the charged pions. The estimations of the cor- 
responding critical field strengths for the charged 
and neutral p mesons are given in Eqs. ([£]) and 



(111. Since these critical strengths are smaller than 
the critical superconducting field ([5| , the supercon- 
ducting condensate should be intrinsically stable, 
at least at the scale of the strong interactions. 

3. The transitions between the unstable and stable 
regions of the p mesons are expected to be smooth 
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crossovers while the onset of the superconductivity 
is expected to be a second order phase transition. 

4. The superconducting state is anisotropic: the elec- 
tric resistance is zero only along the axis of the 
magnetic field. 

5. The superconducting state is inhomogeneous: the 
condensate shares similarity with the Abrikosov 
vortex lattice in the mixed state of a type-II su- 
perconductor. 

6. The pure homogeneous superconducting state is 
not formed. 

7. The onset of the superconductivity of the charged 

mesons leads to emergence of an inhomogeneous 
superfluidity of the neutral p° mesons. The super- 
fluidity is induced by the inhomogeneities of the 
superconducting condensate. 

8. The inhomogeneous superconducting state is real- 
ized as the p-vortex lattice. Locally, the p-vortex 
core expels both superconducting and superfluid 
condensates of the charged and neutral p mesons, 
respectively. The magnetic field takes its maxima 
outside the vortices, while the the strength of the 
superfluid (electrically neutral) field is peaked at 
the vortex centers. However, the unit p-vortex cell 
carries one unit of the quantized magnetic flux of 
the magnetic field while the net p° "flux" is vanish- 
ing. 



9. The spontaneous emergence of the superconduct- 
ing condensate "locks" the rotations of the system 
around the magnetic field axis with a global sub- 
group of the gauge transformations. The inhomo- 
geneities of the condensate break the locked group 
further to the group of discrete rotations of the vor- 
tex lattice. 

10. The Meissner effect (understood in the usual sense) 
cannot be realized in the superconducting QCD 
state due to the Lorenz-invariance of the vacuum. 

Our results also imply that the inhomogeneous 
Ambj0rn-Olesen state [T21 H3] of the vacuum of the elec- 
troweak model is, in fact, an anisotropically supercon- 
ducting state. This state may in principle be realized in 
first moments of the Universe if strong enough magnetic 
fields are created in the primordial era 30 . The su- 
perconducting nature of the Ambj0rn-Olesen state may 
have imprints in the large-scale structure of the magnetic 
fields in the present-day Universe. 
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